Mon. Not. R. Astron. Soc. 000,[Tli](2006) Printed 5 February 2008 (MN KTeX style file v2.2) 



Is Ursa Major II the Progenitor of the Orphan Stream? 
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ABSTRACT 

Prominent in the 'Field of Streams' - the Sloan Digital Sky Survey map of substructure in the 
Galactic halo - is an 'Orphan Stream' without obvious progenitor. In this numerical study, we 
show a possible connection between the newly found dwarf satellite Ursa Major II (UMa II) 
and the Orphan Stream. We provide numerical simulations of the disruption of UMa II that 
match the observational data on the position, distance and morphology of the Orphan Stream. 
We predict the radial velocity of UMa II as —100 kms~^, as well as the existence of strong 
velocity gradients along the Orphan Stream. The velocity dispersion of UMa II is expected 
to be high, though this can be caused both by a high dark matter content or by the presence 
of unbound stars in a disrupted remnant. However, the existence of a gradient in the mean 
radial velocity across UMa II provides a clear-cut distinction between these possibilities. The 
simulations support the idea that some of the anomalous, young halo globular clusters like 
Palomar 1 or Arp 2 or Ruprecht 106 may be physically associated with the Orphan Stream. 

Key words: galaxies: dwarfs — galaxies: individual: UMa II — galaxies: kinematics and 
dynamics — galaxies: evolution — methods: N-body simulations 



1 INTRODUCTION 



Data from the Sloan Digital Sky Survey (SDSS: lYork et al.ll200Q h 
have revealed abundant examples of streams and substruc ture in 
the Milky Way halo. For example, iBelokurov et"aD ( l2006ah used a 
simple colour cut ^ — r < 0.4 to map out the distribution of stars 
in SDSS Data Release 5 (DR5). The "Field of Streams", an RGB- 
composite image composed of magnitude slices of the stellar den- 
sity of these stars, showed the leading arm of the well-known Sagit- 
tarius stream and the Monoceros ring very clearly. Also prominent 
was a new stream, which did not have an identified progen itor, and 
was called the Orphan Stream by IBelokurov et al.l (l2006ah . 

The Orphan Stream was then analysed independently by 
two groups. Grillmair (2006) reported that there was a diminu- 
tive Galactic satellite that lay near the projected path of the new 
stream but that it was "unlikely to be related to it". IBelokurov et al.] 
(l2006bh disagreed, noting that the diminutive satellite lay on the 
same Galactocentric great circle as the Orphan Stream. They ar- 
gued that there was a preponderance of unusual objects along this 
great circle - including the Complex A High Velocity Clouds and 
the young halo globular clusters Ruprecht 106 and Palomar 1 - and 
suggested that some or all may be the remnants of the disruption of 
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a dwarf galaxy. Then JZucker et al.l ( l2006h provided follow-up Sub- 
aru imaging of the diminutive satellite, confirming it as a disrupted 
dwarf galaxy and naming it Ursa Major II (UMa II) after its host 
constellation. 

One possible interpretation of the data is that UMa II is the 
progenitor of the Orphan Stream. Closely related is the possibility 
that both UMa II and the Orphan Stream are remnants from the 
brea k-up of a still larger object, perhaps a tidal dwarf galaxy (see 
e.g.,'Kroupa 1997). In this theoretical study, we strengthen the case 
for such interpretations by providing an orbit for the disruption of 
UMa II so that its tidal tails match the observational data available 
on the Orphan Stream. 



2 OBSERVATIONAL DATA 



The UMa II dwarf galaxy ( IZucker et al.l2006h is located at right as- 
cension a — 132.8° and declination 5 — +63.1°. Its heliocentric 
distance is estimated as Dq = 30 ± 5 kpc, whilst its radial veloc- 
ity is as yet unmeasured. UMa II appears elongated along lines of 
increasing right ascension with an ellipticity of ~ 0.5. Follow-up 
observations of the central parts with Subaru reveal more than one 
density enhancement within the satellite, which supports the fact 
that it may be in the process of tidal disruption. But, no obvious 
tails around the object are discernible in the wider field SDSS data. 
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Figure 1. Orbit of UMa II in the (x, y)-plane (left), (x, zj-plane (middle) and {R = yx^+y^^ 2;)-plane (right). The red solid line is the backwards orbit 
from the present position and the green dashed line is the forward orbit over 1.5 Gyr. The red star with error bar shows the present position of the UMa II dwarf 
galaxy. Black crosses with error bars show the position of the Orphan Stream from Table [T] The positions of some globular clusters which may be associated 
with the stream are marked with blue, open triangles. The distance bracket to the Complex A is marked with black, open squares. 



Table 1. Positions, distance moduli, distances and heliocentric velocities of 
the Orphan Stream From Belokurov et al. (2006b). 



a 5 m — M Dq vq 



162.1° 


-0.5° 


16.5 ±0.7 


20tl 


-35±10kms-i 


158.9° 


8.5° 


16.5 ±0.9 


20l^° kpc 




155.4° 


17.0° 


17.1 ±0.7 


26l^° kpc 




152.3° 


25.0° 


17.5 ±0.8 


32+}^ kpc 




149.4° 


32.0° 


17.5 ±0.9 


32+}^ kpc 


±105 ± lOkms-i 



The total luminosity of UMa II is Mtot,v = -3.8 ± 0.6 mag. This 
translates into a stellar mass of 6 x 10^ M©, applying a conserva- 
tive mass-to-light ratio of 2 which is typical for an old population. 
This is a lower limit for the present-day mass of the remnant object. 

UMa II lies on the same great circle as the Orphan Stream. 
This can be traced for over ^ 50° in upper main sequence and 
turn-off stars in the SDSS data. By constructing a colour-magnitude 
mask ba sed on the ridge- l ine of t he old metal-poor globular clus- 
ter M92. iBelokurov et all (l2006bh showed that the Orphan Stream 
is closer to us at lower declinations than at higher. The distances 
and distance moduli to the Stream at different right ascension and 
declination are listed in Table [T] 

The total magnitude of the Orphan Stream is rrir ~ 9.8. As- 
suming the smallest distance modulus from Table [T] of m — M — 
16.5, this results in Mr ^ —6.7 or 3.5 x 10^ solar luminosities. 
Taking the largest distance modulus of m — M = 17.5, the total 
luminosity of the stream is ^ 8 x 10^ L©. With a mass-to-light 
ratio for an old stellar population of 2, this amounts to a total mass 
in stars in the Orphan Stream of^ 10^ M©. 

IBelokurov et alj (l2006bi) speculated that there might be a con- 
nection between the Orphan Stream and the agglomeration of high 
velocity clouds known as Complex A, which lie on the same great 
circle. Complex A is located between a = 126.7°, 5 = 67.4° 
and a — 134.5°, 5 — 61.7° with a distanc e bracket between 4 
and 15 kpc in heliocentric distance (see e.g., IWakker et al .11 19961 : 
IWakkei]l200lh . The measured radial velocity of this cloud complex 
is in the ranse of -140 to -190kms"\ Of course, the velocities 
of gas clouds may be affected by forces other than gravitational 
ones. 



3 SET-UP 

Our working hypothesis is that the UMa II dwarf galaxy is the 
progenitor of the Orphan Stream. To determine a possible orbit, 
we first perform test-particle integrations in a Milky Way potential 
which consists of a logarithmic halo of the form 

^haio(r) = !|ln(x'±/^-2+d^), (1) 

with q — I^vq — 186 kms~^ and d — 12 kpc. The disc is repre- 
sented by a Miyamoto-Nagai potential: 

^^UR.z) = . (2) 

B? + {h + ^/£^Td^f 

with Md = 10" M0, h = 6.5 kpc and c = 0.26 kpc. Finally, the 
bulge is modelled as a Hemquist potential 

'fbui,e(r) = (3) 

r ± a 

using Mb = 3.4 X 10^° M© and a = 0.7 kpc. The superposi- 
tion of these components gives quite a good representation of the 
Milky Way. The circular speed at the solar radius is 220 km s~^. 
The major advantage is the analytical accessibility of all quantities 
(forces, densities, and so on). 

First, we use trial and error to find a suitable orbit which re- 
produces most of the observational data. We then compute this or- 
bit backwards for 10 Gyr an d insert a live progeni tor. We use the 
particle-mesh code Superbox (iFellhauer et"^l200Qh to perform the 
forward integration until the position of UMa II today is reached. 
We then analyse the location of the tidal tails, adjust the parameters 
from the test-particle simulation and re-run the full N-body model 
to optimise the fit to the observational data. This procedure has to 
be done because the location of the tidal tails differs from that of 
the orbit. 

At outset, we do not distinguish between dark and luminous 
matter and use a one-component model with a Plummer profile. 
Later, we also use a more elaborate two-component model, moti- 
vated by the endpoints of cosmological simulations. It has a Hern- 
quist sphere corresponding to the luminous matter, embedded in a 
dark matter halo which has the Navarro-Frenk- White form. We in- 
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Figure 2. All- sky view of the UMa II orbit and the best matching simu- 
lation. The panels show right ascension versus declination (upper), right 
ascension versus heUocentric distance (middle) and right ascension ver- 
sus heliocentric radial velocity (lower). The grey-scale contours show the 
logarithmic densities of the UMa II tidal tails in the simulation. The solid 
red (dashed green) line is the backward (forward) orbit shown for 1.5 Gyr. 
Black crosses mark the position of the Orphan Stream and the red star marks 
the position of UMa II. Blue triangles show the globular clusters, Arp 2, 
Terzan 7, Ruprecht 106 and Palomar 1. Squares show the position of Com- 
plex A with its distance and velocity brackets. [The starting mass of the 
Plummer model representing UMa IIis4xlO^M0 and its scale-length is 
80 pc] 



vestigate the effects of changing the initial mass, dark matter con- 
tent and scale-lengths of both models. 



4 THE ORBIT 

4.1 Predicted Velocities 



Our best matching orbit puts UMa II at a heliocentri c distance of 
34 kpc , which agrees with the observational datum of IZucker et alj 
(l2006h . We predict the radial velocity and (heliocentric) proper mo- 
tions of UMa II as 



/ila COS S = 



-100 kms"\ 
-0.33 mas yr~^ 



(4) 



IJ.S — —0.51 mas yr 

The resulting orbit is shown in Fig. \T\ It has a perigalacticon of 
^ 18.4 kpc and an apogalacticon of ^ 40.6 kpc. This orbit not 
only connects UMa II with the Orphan Stream but also permits 
Complex A and several globular clusters of the Milky Way to be 
related to it. 

The tidal tails of UMa II do not lie precisely along UMa II's 
orbit. Fig. |2] shows grey-scale contours of the tidal debris in the 
planes of right ascension versus declination, heliocentric distance 
and heliocentric velocity respectively. The positional data on the 
Orphan Stream is nicely matched by the tidal tails. The model 
predicts a strong velocity gradient along the Orphan Stream with 
the radial velocity varying from 200 kms~^ at the southern end 
(a ^ 170°) to -100 kms"^ at the northern end (a ^ 130°). 
The gradient in radial velocity becomes shallower at higher decli- 
nations. 

The Orphan Stream may also have been detected as a density 
enhancement in star count data derived from CADIS or the Calar 
Alto Deep Imaging Survey (iFuchs et al.ll2006h . This idea receives 
some support from Fig.|2l as their 9h field falls on the second wrap 
of the backward orbit. 



4.2 Possibly Associated Objects 

On the basis of intersections of their polar paths, iBelokurov et alJ 
(l2006bh speculated that there may be a connection between the Or- 
phan Stream and a number of anomalous, young halo globular clus- 
ters - in particular Palomar 1, Ruprecht 106, Arp 2 and Terzan 7. 

From Figs. \T\ and |2l we can assess how Belokurov et al.'s 
speculations fare against the simulation. The position and radial 
velocity of Pal 1 is a good match to the forward orb it of UMa II. 
In this context, it is interesting to note that Figure 1 of'Zu cker et al.l 
12006) shows clumps visible in the central parts of UMa II. Pal 1 
looks like one such clump that has already broken off and leads 
UMa II. Arp 2 is also well-matched in position and radial velocity 
of the forward orbit, although it has also been claimed as a possible 
Sagittarius stream member on the basis of distance, kinematics and 
chemical composition (see e.g., Sbordone et al. 2005). The position 
of Rup 106 is a good match to the backwards orbit, but its veloc- 
ity is not (it should lie on the upper rather than the lower wrap in 
Fig. 12]). However, bearing in mind the distance errors to the globular 
clusters, Rup 106 probably cannot be discarded. Ter 7 does seem 
to be ruled out - the right panel of Fig. [T] and the middle panel of 
Fig.|2]show substantial mismatches between its distance and that of 
the forward orbit of UMa II. 

Belokurov et al ] (l2006bh also pointed out the remarkable 
alignment between the Orphan Stream and the Complex A associa- 
tion of High Velocity Clouds (HVCs). Although we do not address 
the origin of Complex A in this paper, we note that the forward 
orbit does pass through the location of Complex A, and even the 
heliocentric velocities are reasonably well-matched (see the lower 
panel of Fig.|2]). If the clouds of Complex A are indeed associated 
with the Orphan Stream, the simulation suggests that they lie more 
than a revolution ahead in orbital phase. 



5 THE MORPHOLOGY OF THE ORPHAN STREAM 

With this orbit in hand, we can deduce some constraints on the 
initial mass of UMa II. The length of the tidal tails is controlled by 
the total initial mass of the satellite (dark and luminous matter are of 
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Figure 3. Left: The length of the tidal tail as a function of the initial mass of the object, shown in the top left of each panel. The red star shows the position of 
UMa II, the black crosses the positions of the Orphan Stream. If the initial mass is of the order 5 x 10^ Mq or more, we should see more than one wrap. If 
the initial mass is of the order 10^ Mq or less, the tail is not long enough to match the Orphan Stream. [The scale-length of the Plummer model representing 
the progenitor of UMa II is 80 pc. The duration of the simulation is 10 Gyr.] Right: The length of the tidal tails as a function of time, shown in the top left of 
each panel. The time has to be on the order of 7.5 Gyr or greater to ensure that the tails are long enough to match the Orphan Stream. [The starting mass of 
UMa II is 4 X lO^M© and its scale-length is 80 pc] 



course not differentiated in our one-component simulations). If this 
mass is < 10^ M©, the resulting tails are too short to be consistent 
with the ^ 50° arc of the Orphan Stream visible in SDSS. On the 
other hand, if the initial mass is ^ 5x10^ M©, further wraps of the 
leading and trailing arms should then be seen in SDSS data. This 
is illustrated in the left panels of Fig. [3] which shows the tidal tails 
produced by the disruption of a sequence of UMa lis of different 
starting masses. Similarly, sequences of the disruption of UMa lis 
for different times, as shown in the right panels of Fig. [3] suggest 
that timescales less than 7.5 Gyr are insufficient to reproduce the 
present-day length of the Orphan Stream. 

Having found lower limits for the progenitor mass and the 
simulation time, we now focus on two particular models. The first 
is a one-component model in which dark and luminous matter are 
not distinguished. It has a Plummer distribution with a mass of 
Mpi = 4 X 10^ Mo and a scale-length of R^i = 80 pc. The ra- 
tionale for these parameters will become clear in Sect. (6] where 
we place further constraints on the progenitor mass by investi- 
gating the remnant. The second is a two-component model with 
the luminous matter represented by a Hernquist sphere with mass 
Mhern = 5x10^ M© and a scale-length of 200 pc. The is embed- 
ded in a Navarro-Frenk- White (NFW) dark matter halo. The NFW 
model has the same scale-length as the luminous matter, together 
with a mass within the cut-off radius (set to be the tidal radius at 
perigalacticon) which is ten times greater than the luminous mass. 
If the NFW mass is made larger or the scale-length smaller, then the 
progenitor becomes much harder to disrupt and does not resemble 
the present-day UMa 11. 

Fig. |4] shows a close-up of the simulation data at the posi- 
tion of the Orphan Stream for the one and two-component models. 



There are a number of morphological features that both the sim- 
ulations reproduce successfully. First, the tidal tails of the models 
have a full- width half-maximum (FWHM) of ~ 2°. This matches 
the FW HM of the Orphan Stream as measured bv lBelokurov et al.l 
(l2006b ). In both simulations, the mass in the Orphan Stream is 
< 10^ M0, in reasonable agreement wit h the stellar mass inferred 
from its luminosity of 8 x lO'^L© (iBelokurov et aLll2006bD . 
There is just one arm visible in the one-component model, and the 
total mass in the Stream is ~ 6 x 10^ M©. This is a closer match 
than the ~ 3 x 10^ M© in stars present in both arms in the Orphan 
Stream for the two-component model. 

Both models reproduce the positional data of the Orphan 
Stream very well and are in good agreement with the measured dis- 
tances. Nevertheless, in the two-component model, a wrap-around 
of the leading arm is present, which gives a better fit to the two low 
declination data points in the middle panels of Fig. |4] Further, the 
velocity data-point at the low declination (or high right ascension) 
end of the Stream shown in the lower panels of Fig. |4] can only be 
reproduced with the presence of a wrapped around leading arm. 



6 THE MORPHOLOGY OF UMA H 

We can sharpen the constraints on the initial mass by requiring that 
the simulations also reproduce the disrupted nature of UMa II itself. 
Fig. [5] shows the results of the disruption of the one-component 
model (left panels) and the two-component model (right panels). 
The three rows show the surface brightness, the logarithmic den- 
sity distribution in right ascension-heliocentric distance and in right 
ascension-heliocentric velocity space, respectively. 

For one-component models, we find that satellites with an 
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Figure 4. A close-up of the Orphan Stream in the one-component model 
(left panels) and the two-component model, showing luminous mass only 
(right panels). Red solid (green dashed) line shows the backward (forward) 
orbit of UM a II. The purple crosses w ith error-bars mark the observational 
results from Belokur ov et ^ l2006b^. From top to bottom the panels show 
the surface brightness in V (mass is converted into luminosity using a mass- 
to-light ratio of 2), the logarithmic density distribution in right ascension- 
heliocentric distance space and the logarithmic density distribution in right 
ascension-radial velocity space. Both models match the positional data of 
the Orphan Stream. Both models also fit the observational distances within 
the errors, but in the two-component model the closest two data-points 
are better matched with the wrap-around of the leading arm, which is not 
present in the one-component model. Also, the velocity measurements are 
only matched if the leading arm is present. 



initial mass > 10^ M© do not become sufficiently dissolved to 
resemble the present-day UMa IL Below this, there is a trade- 
off between starting mass and scale-length. For example, a Plum- 
mer sphere with mass Mpi = 5 x 10^ M© and scale-length 
i^pi = 100 pc gets completely dissolved without a remnant, whilst 
one with — 85 pc gives a remnant which is too massive by two 
orders of magnitude. Reducing the mass to Mpi = 4 x 10^ M© 
and using — 80 pc results in a remnant with similar mass and 
aspect to UMa II. Fig. [6l shows the evolution of the bound mass 
of our one-component model. If the final mass of the remnant is 
~ 6 X 10^ M©, then an object with an initial mass of 10^ M© 
must be in its final stage of dissolution. A robust result is that the 
initial distribution of the satellite cannot be very concentrated, oth- 
erwise there is insufficient mass loss to produce the Orphan Stream. 

For two-component models, it is a challenge to reproduce the 
dissolved nature of the present-day UMa IL As our illustrative ex- 
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Figure 5. A close-up of the UMa II remnant in both models. Left panels 
show again the one-component model, while right panels show the two- 
component model. The top row shows the surface brightness of all parti- 
cles in the one-component model and of the luminous matter only in the 
two-component model. The second row gives the logarithmic density dis- 
tribution in right ascension-distance space while the third row shows the 
distributions in right ascension-radial velocity. The coloured lines are as in 
Fig.a 




t [Gyr] 

Figure 6. The bound mass of UMa II is plotted against time for the one- 
component model. During most of its lifetime, the mass decreases only 
slightly with each perigalacticon passage. However, the last disc shock leads 
to the final disruption of the object and the bound mass drops quickly to 
zero. At this particular instant, the stars of the object become unbound but 
have not yet dispersed from the location of the object. At later epochs, the 
bound mass is zero and the stars disperse into the tails. 
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Figure 7. The peak surface density Emax (black) and the mean heliocen- 
tric distance Dq (green) of the Orphan Stream are plotted as a function 
of right ascension for the one-component (two-component) models in the 
upper (lower) panels. The vertical red lines mark the range of right ascen- 
sion over which the Orphan Stream is detected in SDSS data. UMa II cor- 
responds to the sharp density peak at right ascension a = 132.8°. The 
fading of the Orphan Stream just before it approaches UMa II is caused by 
the decreasing peak surface density and the increasing mean distance. 



ample, we use a Hernquist sphere of 5 x 10^ M© and scale-length 
of 200 pc, embedded in a Navarro-Frenk- White halo with the same 
scale-length and with a mass of 5 x 10^ M© within the tidal radius 
of 400 pc. These parameters are set given the constraint that the 
present-day mass in the Orphan Stream is ^ 10^ Mq. If the mass- 
to-light ratio of the progenitor is ^ 10, this fixes the halo mass, 
whilst the scale-lengths must be in excess of 200 pc to allow for 
enough luminous matter to be stripped off and found in the Orphan 
Stream. Even so, at the endpoint of the simulation, the remnant has 
^ 10^ Mq in stars. This is too large by two orders of magnitude! 

For both the one and two component models, the UMa II rem- 
nant in Fig. \5\ shows a prominent elongation - not along its orbit 
- but along lines of constant declination. The same elongation is 
found in the deepe r, follow-up observ ations with the Subaru tele- 
scope reported by Zuckeretal Comparing the size of the 
observed UMa II of about one degree along constant declination 



Figure 8. Histogram of heliocentric distances of stars in a 0.4° x 0.4° 
field centred on U Ma II, approximately the same size as the panels in 
"Zucker et al.' ^2006). The inset shows the histogram of distances but now 
confined to the very central parts of UMa II (0.1° x 0.1° field). The data 
are taken from the one-component (two-component) model in the upper 
(lower) panels. 



with our models, we conclude that our one-component model fits 
the extension of the real object much better than the two-component 
model. 

Another advantage of the one-component model is that there is 
some substructure in the UMa II remnant, as is visible in the upper 
left panel of Fig. [5] In the simulation, this is caused by tidal shock- 
ing of the remnant at the last few disk passages. The substructure is 
qualitatively similar to internal dumpiness of the UMa II dSph seen 
by Zucker et al. ( 2006). However, to recover the details of this fea- 
ture may well require a more elaborate starting model than a simple 
Plummer sphere. The two-component model shows no substructure 
at the end of the simulation. 

Both simulations not only match all the available positional 
data, but - more strikingly - they also explain why the tails around 
UMa II are faint and undetectable with SDSS. At the positions of 
the Orphan Stream, projection effects enhance the visibility of the 
well-collimated stream, which lies almost along the line of sight. 
By contrast, at UMa II the orbit is almost transverse to the line 
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Figure 9. This shows contours of the V band surface brightness, the velocity dispersion and the mean radial velocity of the remnant in our simulations. The 
first three rows refer to the one-component model after 9 Gyr, 10 Gyr and 11 Gyr. This is a sequence from bound through disrupting to almost completely 
dissolved object. The final row shows the endpoint of the two-component model after 10 Gyr for comparison. (The key to the colour code is given on the 
right-side of the top panels). 
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[kms-1] 


[mas yr~^] 


[mas yr~^] 


[kpc] 


[kpc] 


Miyamoto-Nagai disc, Logarithmic halo 


0.90 


30 


-100 


-0.50 


-0.50 


18.9 


34.5 


Miyamoto-Nagai disc. Logarithmic halo 


1.00 


30 


-115 


-0.40 


-0.50 


18.8 


37.2 


Miyamoto-Nagai disc. Logarithmic halo 


1.11 


30 


-105 


-0.35 


-0.50 


17.8 


36.7 


Dehnen & Binney models 


1.00 


33 


-125 


-0.25 


-0.65 


14.2 


40.8 



Table 2. Parameters for best fit test-particle orbits for different choices of potential. The first column gives the Galactic model. For the logarithmic halos, q is 
the flattening of the equipotentials [see eq. {T)], whereas for the Dehnen & Binney (1998) models, q is the flattening of the isodensity contours [see eq |3]. 
The columns give the best-fit initial conditions of UMa II today so as to join up with the Orphan Streaml; the heliocentric distance, heliocentric radial velocity 
and proper motion in a and 6. The last two columns show the peri- and apogalacticon distances of the orbit. Note that the values for the g = 1 logarithmic 
halo case differ slightly from the values in the main paper because only the test-particle orbit was fit to the data in this Table. 



of sight and there is no enhancement from projection effects. This 
provides a natural explanation as to why an extension of the Or- 
phan Stream is not visible all the way up to the position of UMa II 
in SDSS data. This phenomenon is also illustrated in Fig. |7] which 
shows the peak density in the simulated Orphan Stream dropping, 
and the mean heliocentric distance increasing, as UMa II is ap- 
proached. 



Fig. [8] shows the distribution of heliocentric distances of stars 
in a 0.4° x 0.4° field centred on UMa II remnant, together with an 
inset that records the same information but now confined to the very 
innermost 0.1° x 0.1° field. Although the innermost parts are quite 
confined, the entire object has a significant depth along the line of 
sight of 1 kpc, particularly in the one-component model. This is 
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of the right order of magnitude to cause the b roadening of feature s 
of the colour-magnitude diagram discerned bv lZucker etaPbOOd) . 



7 VELOCITY DISPERSIONS 

We have shown that the initial mass (stars and possible dark matter) 
must exceed 10^ M© to account for the length of the tidal tails and 
the known stellar mass in the Orphan Stream. But any object with a 
total mass > 10^ Mq typically leads to a present-day UMa II which 
is still strongly bound and has a luminous matter contribution at 
least an order of magnitude larger than the observed 6 x 10^ Mq. 

One solution to this dilemma is to postulate that we are ob- 
serving UMa II at a time close to its disintegration. The last disc 
passage led to the almost complete disruption of the remnant ob- 
ject. The stars are now rapidly becoming unbound. They have not 
yet dispersed along the orbit into the tidal tails and we still see them 
in the innermost parts of UMa II in a very confined area (see the in- 
set of Fig. [7]). If so, then the interpretation of kinematic data may 
need special care. 

The first three rows of Fig. [9] all show one-component mod- 
els. We used exactly the same set-up, but started the simulation at 
9 Gyr, 10 Gyr or 11 Gyr ago on the same orbit to make sure all 
models are now seen at the same position on the sky (to exclude 
projection effects). The bound object in the first row has a small 
velocity dispersion, which is even lower than that in the surround- 
ing tails. But the mean line of sight velocity is constant throughout 
the bound object and a gradient is only visible in the tails. This 
changes dramatically in the disrupting model shown in the middle 
panel. We still see an object with a similar total surface brightness, 
but it already shows sub- structure on small scales. The velocity dis- 
persion is inflated by a factor of ten, but the dissolved nature of the 
remnant is already visible in the mean radial velocity. There is a 
strong gradient throughout the object, even though the mean ve- 
locity shows some flocculent structure. In the third row, in which 
the process of disruption is almost complete, the dissolved object 
has a low surface density, which will decrease further in the future 
until it matches that of the tails. The velocity dispersion is again 
low, at much the same value as that of the tails. Looking at the 
mean radial velocity, it is hard to distinguish what remains of the 
object from the tails. For comparison, the final row of Fig. [9] shows 
the same quantities for the two-component model. The final bound 
object has smooth surface brightness contours, a high velocity dis- 
persion because of the dark matter content and no gradient in the 
mean radial velocity. 

Follow-up high precision kinematic observations of this new 
dwarf galaxy could reveal a high velocity dispersion, irrespective 
of the dark matter content. However, the existence of a gradient in 
the mean radial velocity provides a clear-cut distinction between a 
disrupting object and a bound, dark matter dominated object. 



8 THE GALACTIC POTENTIAL 

Hitherto, our Galactic model is built from three fairly simple ana- 
lytic components that could have some deficiencies. Although we 
are using the same standard model as ma ny previous investiga- 
tors (e.g. jHelmill2004l : |johnston et al.ll2005h , it is prudent to exam- 
ine the robustness of our results to changes in the underlying Galac- 
tic potential. Whilst re-running all the N-body simulations would 
be time-consuming, it is straightforward to carry out the initial test- 
particle calculations described in Section 2 for different potentials. 



For example, we can vary the flattening of the halo. If we 
change the halo shape from spherical to moderately prolate or 
oblate, we still are able to fit all the data by slightly altering the 
starting velocities of UMa II. This is illustrated in Table [2l which 
gives the velocities, and the pericentric and apocentric distances 
for test particle calculations. Note that q in the logarithmic halo 
refers to the flattening of the equipotentials - the flattening in the 
density cont ours is typically two or three times greater (see e.g., 
Evan sl ll993h . For moderate changes, a suitable orbit can always be 
found that joins up UMa II with the Orphan Stream, but UMa II's 
predicted velocity and proper motions are then somewhat differ- 
ent. Only if we use strongly prolate or oblate models does the orbit 
of UMa II change so dramatically that we are not able to fit all the 
data on the Orphan Stream at once. However , the recent study of the 
multiple wraps of the Sagittarius' stream bv lFellhauer et al.f ( [200^ 
provides strong evidence that only spherical or close to spherical 
halo shapes are possible for the Milky Way. 

As a further check, we change the type of the Galactic poten- 
tial and use a Dehnen & Binney (1998) model. These potentials 
consist of three exponential discs (thin, thick and gaseous). The 
halo and the bulge are represented by two spheroidal distributions 



ps{R,z) 




Here = i?^ + z^q~'^ and q is the axis ra tio in the density, whils t 
the remaining parameters are chosen as in lFellhauer et al.l ( [20061) . 
Table |2] shows how the intitial conditions or the test particle orbit 
change for the spherical case (^ = 1) for comparison. Again, an 
orbit matching UMa II to the Orphan Stream can be found, and the 
changes in the initial conditions in Table |2] give an indication of 
the likely uncertainties in our predictions caused by changes in the 
Galactic potential. 



9 CONCLUSIONS 

We have carried out N-body simulations to model the evolution and 
disruption of the recently discovered dwarf galaxy UMa II. The 
simulations reproduce the available observational data on the Or- 
phan Stream within their error margins. We conclude that UMa II 
is a likely progenitor of the Orphan Stream. We predict the radial 
velocity of UMa II as —100 kms~^. We also predict a strong ve- 
locity gradient along the Orphan Stream with the radial velocity 
varying from ~ 200 kms~^ at the southern end to ^ —100 kms~^ 
at the northern end. 

From the length of the tails and the mass found in the Orphan 
Stream, we deduce that the initial mass of UMa II is in excess of 
10^ M0. But, an object more massive than 10^ M© cannot be dis- 
solved to produce the present day UMa II, at least on the orbit de- 
rived from the observations. Therefore, the initial mass of UMa II 
has to be of the order a few 10^ M©. To reduce UMa II's mass 
through tidal effects to its present value, the distribution of stars 
and dark matter has to be extended. We carried out a suite of sim- 
ulations of the disruption of UMa II with one-component models, 
which have little dark matter beyond that associated with the stellar 
populations, and two-component models with a mass-to-light ratio 
of ~ 10. There are strengths and weaknesses of both sets of sim- 
ulations. Both reproduce the positions and distances of the Orphan 
Stream, but the two-component models are in bette r agreement with 
the admittedly uncertain kinematic data derived by Belok urov et al.l 
(.2006h) . However, the one-component models can provide a much 
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better match to the disrupted nature of UMa II today. The veloc- 
ity dispersion is not a clean test between these two possibilities, 
as we have shown that objects undergoing disruption can have an 
anomalously high velocity dispersion. However, a clear-cut test is 
provided by the mean radial velocity, which should show no gradi- 
ent for dark matter dominated models, but an obvious gradient for 
disrupting models. 

The orbit that w e have derived supports the idea of 
iBelokurov et"aD ( l2006b) that some of the anomalous, young halo 
globular clusters (particularly Pal 1, Arp 2 and possibly Rup 105) 
may be associated with the Orphan Stream. Intriguingly, the posi- 
tion and velocity of Complex A can also be matched, but only if 
it lies a revolution ahead in orbital phase. The association of these 
objects however makes most sense in the picture in which UMa II, 
the young halo globular clusters and Complex A are all fragments 
of a much larger object like a tidal dwarf galaxy. 
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